THE DERIVE - NEWSLETTER #85 


ISSN 1990-7079 

THE BULLETIN OF THE 



USER GROUP 



Contents: 


1 Letter of the Editor 

2 Editorial - Preview 
Philip Todd 

3 Approximating Circular Arcs with Cubic Splines 

Duncan E. McDougall 
15 Diophantine Polynomials (2) 

Heinrich Ludwig 

22 Cubus Simus - The Snub Cube 
Guido Herweyers 

35 Statistics with TI-Nspire (1) 

43 Our USER Forum 

45 TIME 2012 - Lectures and Workshops 


March 2012 








D-N-L#85 


Information 


D-N-L#85 


Neu: ACDCA-Mitgliedschaft 

Das ACDCA (Austrian Center for Didactics of Computer Algebra) wurde Ende der 80er 
Jahre gegriindet und ist seither eine der treibenden Krafte fur den Einsatz von CAS in der 
Mathematikausbildung. Dieser Verein wurde kurzlich umstrukturiert und bietet nun alien 
Interessierten die Mitgliedschaft und damit die Gelegenheit zur Mitwirkung an. Bitte in- 
formieren Sie sich auf der Homepage des ACDCA iiber die Aufgaben und Ziele des 
ACDCA. 

Info und Beitrittserklarung finden Sie auf 

http://rfdz.ph-noe.ac.at/acdca.html 

Hier steht auch die verbesserte und erganzte Version der „Dynamischen Systeme“ zum 
Download bereit: 

http://rfdz.ph-noe.ac.at/acdca/materialien.html 

Latest News: 

You are invited to download 

Dynamic Systems on Various Platforms - An Excursion from 
Environment & Tourism to Strange Attractors! (113 pages) 

http://rfdz.ph-noe.ac.at/acdca/materialien.html 


Prof. Paditz from HTW Dresden informed that the complete proceedings of two MEC21- 
Conferences (Dresden 2009 with 147 contributions and Grahamstown 2011 with 73 con- 
tributions) are available as download at 

http://nbn-resolving.de/um:nbn:de:bsz:14-qucosa-79236 

http://nbn-resolving.de/um:nbn:de:bsz:14-qucosa-79243 


STATISTICS with Tl-Nspire 

I am very happy that Guido Herweyers - our Flemish Statistics expert - permitted publi- 
cation of his 50 pages paper. I tried to translate from Flemish (Dutch) into English (my 
wife is from Dutch origin) and you will find his paper divided in some parts in this and in 
the next DNLs. TI-Nspire has some very nice features especially for statistics. 


Attend TIME2012 in Tartu, Estonia 
July 10-14, 2012 
http://time2012.ut.ee/ 


List of lectures and workshops, which have been accepted until now is on pages 
45 and 46. Don’t miss the opportunity to meet CAS-experts from all continents 
and register! 
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Liebe bUG-Mitglieder, 

Leider bin ich mit diesem bNL etwas zu 
spat dran weil ich die Arbeit an einigen 
Artikeln unterschatzt habe und die Vorbe- 
reitung der TIME 2012 immer wieder Zeit 
in Anspruch genommen hat. 

Im letzten bNL wurde das TT-Nspire- 
Programm TRIGO vorgestellt. Herr van 
Lantschoot hat eine verbesserte Version 
geschickt, herzlichen bank dafiir. 

Es freut mich besonders, dass Philip 
Todd einen Artikel beigetragen hat, in dem 
er das Zusammenwirken seines schonen 
Programms Geometry Expressions mit 7T- 
NspireCAS demonstriert. bass dies auch 
mit DERIVE { u.a. CAS) moglich ist, habe ich 
zu zeigen versucht. Philip wird auch an der 
TIME 2012 mitwirken. 

Ludwig Heinrichs Cubus-Simus-Aufsatz 
hat auch fur mich eine Eigendynamik be- 
kommen. Viel Information war im Internet 
zu finden. Einige Challenges" fur Sie sind 
noch off en! 

Ich kann mir denken, dass bietmar 
Oertel schon ungeduldig auf die Veroffent- 
lichung seiner umfangreichen Untersuchun- 
gen im Bereich der Primzahlen und der Tay- 
lorreihen wartet. Seine Unterlagen umfas- 
sen jeweils ca 60-70 Seiten und eine sinn- 
volle Zusammenfassung braucht seine Zeit, 
aber ich bin schon dabei. 

ber im letzten bNL angekundigte Vor- 
schlag von Klaus Korner muss auf den 
bNL#86 warten wie auch interessante 
Informationen fur mu-Math-Kenner. 

Von der Webseite des ACbCA konnen 
Sie eine verbesserte und erweiterte Versi- 
on meines Aufsatzes zu den Dynamischen 
Systemen (ca 110 Seiten) beziehen. Beach- 
ten Sie bitte unsere Einladung auf der In- 
fo-Seite. 

SchlieGlich finden Sie die Liste der bis- 
her angenommenen Vortrage fur die TIME 
2012, in deren Rahmen auch ein bUG- 
Treffen stattfinden wird. 


bear bUG Members, 

It is really a pity that I am a little bit 
late with this bNL. Actually I underesti- 
mated work on some contributions and 
preparation of TIME 2012 needed some 
time, too 

The TI-Nspire-program TRIGO was 
presented in the last bNL. Erik van 
Lantschoot, the author of TRIGO sent an 
improved version, many thanks for this. 

I am very happy that Philip Todd pro- 
vided an article demonstrating cooperation 
of his excellent Geometry Expressions pro- 
gram with TI-NspireCAS. I try to show 
that this is also possible with DERIVE (and 
other CAS). Philip will present a lecture and 
a workshop at TIME 2012. 

Ludwig Heinrich's Cubus Simus contri- 
bution got an interesting self dynamics for 
me. Much information could be found in the 
web. Some “challenges" are still remaining 
for you. 

I am quite sure that bietmar Oertel is 
is anxious for publication of his extensive 
investigations in the field of prime numbers 
and Taylor series. His papers comprise 
about 60-70 pages each, so preparing a 
meaningful summary needs time, but I am 
busy ... 

I announced a proposal from Klaus 
Korner in the last bNL. It must wait for 
bNL#86 - together with interesting infor- 
mation for my-Math-experts. 

You are invited to download a paper on 
Dynamic Systems on Various Platforms 
(~100 pages) from the ACbCA-website. Its 
URL is given on the information page. 

Finally you can find the list of the ac- 
cepted lectures and workshops for TIME 
2012 where we will have a bUG-Meeting, 
too. 

Viele GruBe, Best regards 



Download all DTVL-Derive- and Tl-files from 

http : / / www . austromath . at/ dug/ 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 7/ User Group. 
It is published at least four times a year 
with a content of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE, 
TI- CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 


Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wiirmla 
Austria 

Phone: ++43-(0)660 3136365 

e-mail: nojo.boehm@pgv.at 


I II 

Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS- 77 Newsletter will be. 


Next issue: June 2012 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 
Wonderful World of Pedal Curves, J. Bohm, AUT 
Tools for 3D-Problems, P. Luke-Rosendahl, GER 
Hill-Encription, J. Bohm, AUT 

Simulating a Graphing Calculator in DERIVE, J. Bohm, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Overcoming Branch & Bound by Simulation, J. Bohm, AUT 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - Interesting features, T. Koller & J. Bohm, AUT 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what's new? J. Bohm, AUT 

Truth Tables on the TI, M. R. Phillips, USA 

Where oh Where is It? (GPS with CAS), C. & P. Leinbach, USA 

Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZK 

Tutorials for the NSpireCAS, G. Herweyers, BEL 

Some Projects with Students, R. Schroder, GER 

Structures in the Set of Prime Numbers, D. Oertel, GER 

Dirac Algebra, Clifford Algebra, D. R. Lunsford, USA 

Treating Differential Equations (M. Beaudin, G. Piccard, Ch. Trottier), CAN 

A New Approach to Taylor Series, D. Oertel, GER 

Cesar Multiplication, G. Schodl, AUT 

Henon & Co; Find your very own Strange Attractor, J. Bohm, AUT 
Rational Hooks, J. Lechner, AUT 

Simulation of Dynamic Systems with various Tools, J. Bohm, AUT 

Using DERIVE to simulate the basic steps of the Q. Shor’s algorithm, N. Urrego, COL 

and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
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I am very happy that Philip kept his promise to send a contribution for the 
DNL. I appreciate his Geometry Expressions as an extremely powerful geome- 
try tool which has a strong built-in CAS. GE can cooperate with all main Com- 
puter Algebra Systems. Come and attend Philip’s lecture and workshop at 
TIME 2012. 

Josef 

Approximating Circular Arcs with Cubic Splines 

Philip Todd, Saltire Software 
philt(a), saltire. com 

We use Geometry Expressions and the TI nSpire CAS to find the best 
cubic spline approximation to a quarter circle. 


1. Introduction 

A problem came from a user of Geometry Expressions recently, which nicely illustrates the power of a 
symbolic geometry system used in conjunction with CAS. The user wanted to find the control points 
of a Bezier spline which approximated a circular arc of a given angle. Why would he want to do such 
a thing? Some software and some machine tools have Bezier spline curves as primitives, but do not 
have circular arcs. If you want a circle, you need to piece an approximation together out of splines. 



Figure 1: A cubic spline constructed in Geometry Expressions. The spline ’s para- 
metric equation is computed by the software 

A Bezier spline can be readily modeled in Geometry Expressions [1] using the geometric version of de 
Casteljau’s algorithm [2]. We start with control polygon ABCD. These points are constrained to have 
coordinates (x 0 ,yo), (xi,yi), (x 2 ,y 2 ), (x 3 ,y 3 ). Point E is constrained to be proportion t along AB. Point F 
is constrained to be proportion t along BC, and point G is constrained to be proportion t along CD. 
Points H and I are now added and constrained to be proportion t along segments EF and FG. Finally, 
point J is added and constrained to be proportion t along segment HI. The spline curve is the locus of J 
as t runs from 0 to 1 . 

An important property of the spline curve is that it is tangential to the control polygon at its ends. That 
is at t=0, the spline curve is tangential to AB, and at t=l it is tangential to CD. 
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This property can be readily verified using the TI-N spire CAS. Geometry Expressions will automati- 
cally generate the parametric equation of the curve (figure 1), and this can be copied (using Copy As / 
TI Nspire) into a TI-Nspire CAS worksheet (figure 2). 


[x=x[o]+x[o] • /• -3+x[o] • f 2 • 3+x[o] • A - l+x[ 1 ] • /• 3+x[ 1 ] ■ f 2 • -6+x[ 1 ] • A 3+x[2] • / 2 • 3+x[2] ■ f 3 • -3+x[3] * 
[x=/ 3 '(x[3]-3-x[2]+3 x[l]-x[o])+/ 2 -(3 x[2]-6 x[l]+3'x[o])+r(3 x[l]-3'x[o])+x[o] jy=f 3 -(y[3]-3^ 

— (r 3 (x[3]-3-x[2]+3-x[l]-x[o])+r 2 -(3x[2]-6-x[l]+3x[o])+/(3x[l]-3-x[o])+x[o])|r=0 
dt 

3'x[l]-3'x[o] 

^(' 3 '(>{3]-3>{2]+3'>{l]-^0])+/ 2 '(3'>{2]-6.y[l]+3'>{0])+('(3.y[l]-3'>{0])+>{0])|/=0 

3>[l]-3^[0] 

— (r 3 (x[3]-3x[2]+3'x[l]-x[o])+r 2 '(3x[2]-6x[l]+3x[o])+/(3x[l]-3x[o])+x[o])|r=l 
dt 

3x[3]-3x[2] 

3 • (y[3]-3 >[2]+3 1 ]->{o])+/ 2 • (3 2]-6 ,y[l ]+3 >(o])+^ (3 1 ]~3 >{o])+>{o])|/= 1 

3'j{3]-3-j{2] 


5/99 

Figure 2: differentiating the spline curve at t=0 and t=l shows 
that it is tangential to the control polygon 


The derivatives of the x coordinate and y coordinate with respect to t evaluated at t = 0 show that the 
tangent aligns itself with the segment AB. Again, the derivatives evaluated at t=l show that this tan- 
gent aligns itself with the segment CD. 


# 1 : 


2 3 2 3 2 3 

x = x + x • t • (-3) + x *t *3 + x *t .(-1) + x »t«3 + x -t • (-6) +x «t *3 + x -t *3+x *t *(-3) + 
00 00 11 122 


3 2 3 2 3 2 

x -t , y = y + y -t.(-3) + y «t «3 + y -t • (-1) + y «t*3 + y *t -(-6) +y -t -3 + y -3+ 
3 00 00 11 12 


3 3 

y -t .(-3) + y -t 
2 3 


d 2 3 2 3 2 3 

#2: — x + x • t* (-3) + x -t -3 + x -t • (-1) + x -t-3 + x -t .(-6) +x -t -3 + x -t -3 + x -t .(-3) + 

dt L 0 0 0 0 11 1 2 2 


3 2 3 2 3 2 3 

x -t , y + y • t* (-3) + y -t -3 + y -t .(-1) + y -t-3 + y -t • (-6) + y -t -3 + y -t -3 + y -t *(-3) + 
300 00 11 122 


#3: 


3' 

y -t 

3 

2 2 2 2 2 2 
3 • t .x + (6* t - 9* t ) • x + (9* t - 12- 1 + 3) • x - (3-t - 6-t + 3)-x , 3-t .y + (6-t - 9-t ).y + 

3 2 1 0 3 2 


#4: 


# 5 : 


2 2 
(9* t - 12- 1 + 3).y - (3-t - 6* t + 3).y 

1 0 . 

[3-x - 3-x , 3*y - 3 • y 1 

L 1 0 1 oj 



DERIVE screen containing expressions imported from GE and further calculations. 
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2. Center Point on Circle 


In this article, we’ll use a spline to approximate a quarter circle, we’ll refer to general arcs briefly in 
the conclusion. 

We model a quarter circle of unit radius in the first quadrant. Point A is (0,1) and point D is (1,0). AB 
is constrained to be perpendicular to the y axis (hence the spline will be horizontal at A). CD is con- 
strained to be perpendicular to the x-axis, and hence the spline is vertical at D. 

It also seems reasonable that we should require our spline to be symmetric about x=y. We enforce this 
by constraining AB and BC to be the same length: k. 



Figure 3: Spline curve with control polygon constrained to guarantee tangency with the unit 
circle at A and D and symmetry about y=x (proportional constraints are hidden for clarity). 


We would like to find a value for k which provides a close approximation to a quarter circle. One 
approach would be to require the center point of the spline curve to lie on the circle. The center point 
is the point with parametric location on the curve Vi. This can be easily created in Geometry Expres- 
sions by constraining a point’s parametric location on the curve to be Vi (figure 4). 

Geometry Expressions can supply the distance of this point from the origin. It hardly needs a CAS to 
find a value k such that this distance is 1. However, as we have one available, we’ll use it. The result- 
ing value of k can be copied from the TI-Nspire CAS back into geometry Expressions. Comparison 
with a unit circle (figure 5) shows that the fit is quite good. How good we will see soon. 



<rm> * *Un$aved ▼ 310 


(4+A: 3) 2 2 1 
8 


( Z k+i) {2 
8 

solve 

f (3 M iE 1; 1 

8 j 

l 


1 AMD. 552285 

► Decimal 


Figure 4: Distance of the center of the spline from the origin is measured in Geometry 
Expressions. The CAS is used to find the value of k which makes this distance 1. 
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(4 + k-3)-2 -1 


#6: 


#7 : SOLUTIONS 


1/2 

(4 + k-3)-2 -1 


= 1, k 



' 4-^2 4 ' 

/ 

3 3 . 


#8: 


[0.5522847498] 



Figure 5: Values of k which make point L lie on the unit circle are copied from the 
CAS into Geometry Expressions allowing a visual inspection of the solution. 

2.1 Error Analysis 

To see how far away from the circle our curve lies, we change the location of L from 14 to a generic 
parametric location x (we’re using x for convenience graphing in the TI-Nspire, but at the risk of con- 
fusion . . .x is not the x coordinate of L, but its parametric location on the curve). Geometry Expres- 
sions generates an expression for the distance of this generic point from the origin (figure below). 
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il+l-408+2 2 •288/ x S +l-204+2 2 • 144/x 3 +\34+2 2 '-24/ x 2 +ll36+2 2 -96/x 6 +\442+2 2 -312J'X 4 J 

J-8 (l2-JT-17) x 6 +24 (l2 ^-17)'X 5 -26 (l2 l /i’-17) x 4 +12 (l2 j2’-17) x 3 -2 (l2 ,/2"-17) x 2 +1 


(J-8-(l2 JT-17) x®+24'(l2 y/T-17) x 5 -26-(l2 ,/2"-17)'X 4 +12 (l2 (/2"-17)'X 3 -2 (l2 ,/T-17) x 2 +1 / 
-[8[\2j2-17)-x 6 -24(l2j2-n}x 5 +26[l2j2-17)-x i -12{l2j2-17)x 3 +2{l2j2-\7)-x 2 - 


— -(8 (l2 i2"-17) x 6 -24 (l2 i2"-17) x 5 +26 (l2 J2"-17) x 4 -12'(l2 j2’-17) x 3 +2 (l2'i2"-17) x 2 -l) 

-4-(l2 ,/2"-17) x (l2 x 4 -30 x 3 +26'X 2 -9 x+l) 


solve(-4-(l2 ,/T-17) x (l2 x 4 -30 x 3 +26 x 2 -9 x+l)=0,x) 


-( x /i’-3) 1 

x=0 or x=— or x=— or x= 

6 2 


Sor x=l 


1/4 


Figure 6: The distance from the origin of point at parametric location x on the spline 
curve is copied from Geometry Expressions into the Tl-nSpire CAS. The expression is 

squared , differentiated and solved. 


To find the values of x for which this distance is maximal or minimal, we first copy the expression 
into our CAS. Optimizing the distance is the same as optimizing its square, so our first step inside the 
CAS is to square the expression, thus getting rid of the square root, and making it more likely that the 
derivative will fit on the page. Next, differentiate, then solve (figure 6). 


1/2 5 1/2 3 1/2 2 

d := (1 + (-408 + 2 -288) -t + (-204 + 2 -144). t + (34 + 2 .(-24))-t + (136 + 

1/2 6 1/2 4 1/2 

2 • (-96)) • t + (442 + 2 *(-312))-t ) 


SOLUTIONS 


fd 2 

— d = 0, t 


0 , 1 , 



1 

2 


V3 


We observe that we have extrema at 0, 1/2 (guaranteed by symmetry), 1 and at two other values. By 
construction, the curve lies on the unit circle at 0, VS and 1, leaving two points of interest, which are 
symmetric about x = VS and should, by the symmetry of the figure, lie at the same distance from the 
origin. This then is the quantity we are interested in. 


We simply copy it from the CAS and paste into the Geometry Expressions , replacing x as the paramet- 
ric location of point L (figure 7). 
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A 



Figure 7: Maximal distance of point L is obtained by pasting the solution forx into the parametric 
location constraint for point L. Numeric display shows the error to be less that 3/1 0,000 



Figure 8: Error graphed as a function of parametric location on the curve. 

An error function may be obtained by subtracting 1 (the radius of the circle) from the distance function 
measured in figure 6. This is displayed in figure 8. We observe that this spline curve always errs on 
the side of being too large. For some applications this may be a benefit, however, we might ask our- 
selves whether it would be possible to improve our estimate if we allowed the spline to err both inside 
and outside the circle. We’ll look at one such curve in the next section. 


3. Equal Area 

One way of ensuring that the spline errs both on the inside and on the outside of the unit circle is to 
insist that the area under the spline is the same as the area under the quarter circle. To measure the 
area under the spline in Geometry Expressions , you need to create first an arc LM on the spline, then 
draw line segments LK and KM (where K is the origin). You can then select the segments and the arc 
and create a polygon. As a final step, you will make sure L and M lie at the start and end of the curve 
by constraining their parametric locations to 0 and 1 . 
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20 5 2 20 5 2 


solve 

f 1 . *' 3 . ,2 - 3 _w 

, -(7-3 (5 ti - 22 ) -6) , J-3'(5 w - 22)+6 


\2 5 20 4 J 

k= or k= 

3 3 

1 , -U- 3 - (5 71 - 22 ) - 6 ) , d - 3 - (5-7t— 22 ) +6^ . , 

\k=-* or k=— ►Decimal 

l 3 3 j 

k= 0.551778 or *=3.44822 




3/99 


Figure 9: The area under the spline is measured in Geometry Expressions by creating an arc LM on 
the spline, connecting L and M to the origin and creating a curve-sided polygon KLM. Finally, L and 
M are constrained to lie at parametric locations 1 and 0 on the curve. In our CAS, we can find the 
value of k which makes the area under the spline equal to the area under a quarter unit circle. 


Copying the area value into the CAS we can equate it with 7i/4 (the area under the quarter unit circle), 
and solve for k. We have two solutions. The first is numerically close to the value of k which puts the 
center point on the unit circle. The second is too far removed to be reasonable. Figure 10 shows what 
they look like: 



Figure 10: two solutions yielding the same area under the spline. Only one is a good 
approximation for the quarter circle! 
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We can paste this value back into Geometry Expressions for k, and let the software derive the distance 
from the origin to the point on the curve at parametric location x (figure 1 1 . This can in turn be copied 
into TI-Nspire for graphing (figure 12). 


6 -^- 3 - (- 22 + 5 - tt) 



Figure 11: distance to the origin for point at parametric location x on the equiareal approximation 

Analysis of the graph shows that while we have succeeded in reducing the extent to which the ap- 
proximate curve overshoots the circle, we have done so at a cost of undershooting it by almost as 
much at the center point. 



Figure 12: Radius for the equiareal approximation is copied into TI-Nspire from Ge- 
ometry Expressions, and plotted next to the area for the central-point approximation. 

However, this suggests an approach for finding a better solution. Can we find a value of k such that 
the error function undershoots by the same amount that it overshoots? 
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4. Minimizing the maximum error 

Using Geometry Expressions , we can generate a distance function for a point at parametric location x, 
copy this into TI-Nspire CAS and differentiate to find the extrema (figure 13,14). 



Figure 13: Distance from the origin of the generic point on the spline with parameter k 


(l+(’24+A:72+A: 2 '-54)-x 5 +(-6+A:6+/r 2 -9)x 2 +(4+A:12+^ 2 '-36)-x 3 +(8+A:-24+A: 2 -18)-x 6 +(l8+A:--66+^ 2 63)'X 4 ) 


2 (9-/: 2 -12'A:+4) x 6 -6 (9^ 2 -12-A:+4)-x 5 +3 (21 J: 2 -22'A:+6) x 4 -4.(9.^ 2 -3.A:-l)-x 3 +3 (3^ 2 +2 ^-2 )-x 2 +1 


(]2-(9^ 2 -12 ^+4) x 6 -6'(9J: 2 -12 '^+4)-x 5 +3 (21^ 2 -22 '^+6) x 4 -4'(9 J: 2 -3 ^-l) x 3 +3-(3^ 2 +2 A:-2)'X 2 +l f 

2(9/: 2 -12-/t+4)x 6 -6(9^ 2 -12A:+4)'X 5 +3(2lJ: 2 -22A-+6)x 4 -4(9-/r 2 -3'A:-l)-x 3 +3-(3^ 2 +2J:-2)-x 2 +l 


— (2-(9^ 2 -12^+4) x 6 -6'(9 '^ 2 -12 A:+4)-x 5 +3'(2r/r 2 -22-A:+6) x 4 -4'(9 ^ 2 -3^-l)-x 3 +3'(3 /t 2 +2 A:-2)'X 2 +l) 


dx 


12 (9-^ 2 -12 ^:+4)-x 5 -30'(9J: 2 -12 A:+4)-x 4 +12 (21 ^: 2 -22 A:+6) x 3 -12-(9J: 2 -3 ^-l) x 2 +6-(3^ 2 +2J:-2)'; 


solve(l2'(9-/r 2 -12A:+4)x 5 -30-(9-/: 2 -12^+4)'X 4 +12(21'/r 2 -22^+6)-x 3 -12-(9^ 2 -3^-l)x 2 +6(3J: 2 +2^-2)-x=0y 


(< 

l-3 * 2 -20 *+12 |3 *-2|-9 * 2 +12 *-4) 

1 

J-3|'* 2 -20-A:+12 • |3- ,A:-2|h-9- /l 2 -12- A+4 

2-(9-* 2 -12*+4) 

■ or x= 

2'(9-* 2 -12-/c+4) 


or x=0 or x=— or x=l 
2 




Figure 14: Extrema for the curve of figure 12 calculated in the Tl-nSpire CAS 


The extrema lie at x = l A and at two more complicated points. We can copy and paste the complicated 
points into Geometry Expressions and let it worry about the algebra. The distance value (figure 15) is 
surprisingly simple. 
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Figure 15: Extreme values for the distance from the origin on the spline curve with parameter k 

In the TI-Nspire CAS we define a variable errorl to be the distance at the midpoint subtracted from 1. 
We define error2 to be the distance at the more complicated extremum minus 1 . The CAS will find, 
numerically, a value of k which equates these errors (figure 16). 


1 

. , (4+£-3)-2 2 '1 
Define error 1=1 

1 8 

Done 


1 



Define error2={~ 1+fc 3) ■ (s+fr -24+ ft? ■ 12+ fr 3, 8+ A 4, 3 ) ■ 

(-2+*'3)" 2 -l 

Done 


so 1 ve ( e rro rl=e rro r2 } ft) 

k= 0.551915 



1 



3/3 


Figure 16: Numerical solution for k, equating the absolute values of the two errors 

A plot of the error functions (figure 17) reveals that the minimax solution reduces the maximal error 
by over 25% compared with either of the other approximations. 
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Figure 1 7: Error curve for the minimax solution lies between the error curves for the 
central point solution and the equiareal solution. 


5. Conclusion 

We have looked at modeling a quarter circle. It is also useful to be able to approximate a general por- 
tion of a circular arc. The development of spline approximations for general arcs yields more compli- 
cated formulas, but you can use a similar process, and the CAS can take care of the details. You do 
need to be careful with the minimax solution as the location of the extremum is complicated enough 
that computing its value can cause problems both for Geometry Expressions and the TI-Nspire. 

A straightforward, but interesting analysis would be to investigate the improvement in accuracy to be 
obtained by approximating a quarter circle by two 45 degree spline approximations. 

Another question would be this: let’s say we are approximating a full circle by splines, each of which 
approximates one quarter of the circle. It is a little artificial to insist that the spline meets the circle at 
the axes. Why not allow the spline to start and finish at a radius slightly below that of the circle being 
approximated. In this way we could essentially translate the error curve for the center point approxi- 
mation vertically down by half its amplitude and achieve a 50% reduction in error. 

The huge benefit of having such powerful software tools as Geometry Expressions and the TI-Nspire 
CAS working together under your direction is that they allow you to explore many more interesting 
variants of a problem than you might be inclined to undertake working only by hand. 

The beauty of Geometry Expressions is that it assists in the mathematical modeling phase of problem 
formulation, allows you to copy the mathematics into your CAS, and then assists in interpreting the 
results, by displaying the solution. In this way degenerate or nonsensical solutions can be immediately 
differentiated visually from sensible solutions. Such was the case in figure 10, where a picture imme- 
diately showed us that one solution involved a spline with a loop. 


References 

[1] www.geometryexpressions.com 

[2] Farin, Gerald & Hansford, Dianne (2000). The Essentials of CAGD.: A K Peters, Ltd. ISBN 1- 
56881-123-3 


p 14 


Ph. Todd: Approximating Circular arcs with Cubic Splines 


D-N-L#85 


Here are some screenshots from my Geometry Expressions - DERIVE collaboration: 


2 3 2 3 2 3 3 

x + x • t • (-3) + x -t -3 + x -t -(-I) + x -t-3 + x -t • (-6) + x -t -3 + x -t -3 + x -t -(-3) + x -t , y 
.00 0 0 11 12 2 30 

2 3 2 3 2 3 3 

+ y -t-(-3) + y -t -3 + y -t *(-l) + y -t-3 + y -t -(-6) + y -t -3 + y -t -3 + y -t -(-3) + y -t 
0 00 11 122 3. 

[ 2 3 2 3 2 3 3 

[o + 0 • t • (-3) + 0 • t -3 + 0 • t .(-1) + k • t -3 + k • t .(-6) + k-t -3 + 1-t -3 + 1-t .(-3) + 1-1,1 + l-t-(-3) 

2 3 2 3 2 3 31 

+ 1-t -3 + 1-t -(-1) + 1-t -3 + 1-t -(-6) + 1-t -3 + k-t -3 + k-t -(-3) + 0-t J 

[32 32 23321 

L k - (3 • t - 6 • t + 3-t) - 2 • t + 3 • t , k • (3 • t - 3-t ) + 2-t - 3-t + lj 

1/2 

(4 + k -3) -2 -1 

error 1 := 1 

8 

2 3 41/2 

(-1 + k • 3) • (8 + - 24-k + 12 • k + 8-k + 3-k ) 

error2 := - 1 

2 

(3-k - 2) 

NSOLVE(er rorl = error2, k) 
k = 0.5519150238 

2 2 

x + y = 1 

[ 3 2 3 2 2 3 3 

Lo . 5519150238- (3 • t - 6-t + 3-t) - 2-t + 3-t , 0 . 5519150238- (3 -t - 3-t ) + 2-t - 


3-t + 1 



The blue curve changes when moving the slider, the black line is the circle line and the spline is plot- 
ted in red color. 

I don’t show plots of the error curves, it is too obvious that this is no problem for DERIVE (and other 
CAS, of course). 
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Diophantine Polynomials 1 ^ (Part 2) 

Duncan E. McDougall, duncanemcdougall@hotmail.com 

( 1 ) 


In the pursuit of families of polynomials which are factorable in the set of integers and whose de- 
rivatives are factorable in the set of rational numbers, it is amazing how much easier the process is 
when we incorporate an old trick from the remainder theorem. I refer to taking the sum of all the nu- 
merical coefficients. If this sum is zero, then the number 1 is an automatic root; if not, it is not. In 
essence, we plan to have one as a root by making sure that the selection of numerical coefficients add 
up to zero in any of the cubic or quartic forms. In terms of the quintic, it does simplify this form to 
some extent, but Galois Theory promises that we will have irrational numbers in exact form. 

To begin, examine all three degrees in various forms with one as an automatic root and obtain new 
families of curves. 


Cubic Polynomials: Form y - ax’ - bx 2 + (b - a)x 


Quickly, a - b + (b - a) is zero and so 1 is root (remainder Theorem). Using synthetic division 


gives 

1 a -b b-a or 

a a-b 
a a-b 0 


(x - 1) ■ (ax 2 + (a- b)x) 

- x • (x - 1) • (ax + (a- b )) 


For integral roots, “b” needs to be a multiple of “a” and this will guarantee that the original poly- 
nomial is factorable. Now take the derivative and look for rational roots; this gives 

y f = 3 ax 2 - 2 bx + (b-a) 

0 = 3 ax 2 -2b-x + (b-a) 


( 2 ) 

Now test the discriminant B 2 - 4 AC to determine values for which the polynomial is factorable. 
Here B = - 2b , A = 3 a, and C = b-a and this gives 

B 2 - 4 AC = (-2 b) 2 - 4(3 a)(b - a) Here 4 is a perfect square, but b 2 - 3ab + 3 a 2 is a 

= 4 b 2 -12 ab + 12a 2 perfect square for a limited number of values. 

= 4(b 2 -3ab + 3a 2 ) 

Fixing values of a in order to obtain values of b to make the expression b 2 - 3 ab + 3 a 2 , a perfect 
square, the following table is produced: 


[1] Diophantine Polynomials (Part 1) appeared in DNL#81. 
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TABLE I 



Value of a 

Resulting Polynomial 

Corresponding Values of b 

1 

b 2 -3b + 3 

1,2 

2 

b 2 - 66 + 12 

2,4 

3 

b 2 - 9b + 21 

-2, 3, 6, 1 1 

4 

b 2 - 126 + 48 

4,8 

5 

b 2 - 156 + 75 

-11,2,5, 10, 13,26 

6 

b 2 -186 + 108 

-4, 6, 12, 22 

7 

b 2 -216 + 147 

-26, -1,7, 14, 22,47 

8 

b 2 - 246 + 192 

1, 8, 11, 13, 16, 23 

9 

b 2 - 276 + 243 

-47,-6, 9, 18, 33 

10 

b 2 - 30b + 300 

-22, 4, 10, 20, 26 


In the interest of discovering single roots only, those values of b which are multiples of a are 
eliminated from consideration because we get double roots. Thus no values of b work for a = 1 or 
a = 2 , but we do get values for b when a = 3 namely b = -2 and b = 11. There are obviously more 
values for b but do not consider them in order to keep the factors relatively small and easy to calculate. 
The above values are easy to spot on the TI graphing calculator by entering y = x 2 - 3x + 3 , and pressing 
2 nd Function Table. 

Simply look for those values of y which are perfect squares. 


( 3 ) 

An example of a cubic can be taken from TABLE I where a = 3 and b =11. This gives 

y - 3x 3 -llx 2 + 8x 
y = x-(3x 2 - llx + 8) 

0 = x • (3x 2 - 1 lx + 8) 

0 = x • (3x 2 - 3x - 8x + 8) 

0 = x • (3x • (x - 1) - 8 • (x - 1)) 

0 = x • (x - 1) • (3x - 8) 


where x = 0, x = 1 and x = 


8 

3 


and for the first derivative see next page. 
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y = 3x 3 - 1 lx 2 + 8x 
y' = 9x 2 - 22x + 8 
y ' = 9x 2 - 1 8x - 4x + 8 
0 = 9jc • (jc — 2) — 4 • (jc — 2) 
0 = (9x - 4) • (x - 2) 


where x = — , x = 2 * 
9 


* This gives rise to a new family of cubics whose original 
form and derivative form are both factorable. 


Polynomial 

x-(x-l)*(3x-8) 

(x - 1) • (x - 2) • (3x -11) 
(x - 2) • (x - 3) • (3x - 14) 
(x - 3) • (x - 4) • (3x - 17) 


Derivative 
(x - 2) • (9x - 4) 
(x-3)-(9x-13) 
(x - 4) • (9x - 22) 
(x-5)-(9x-31) 


)uartic Form 


Moving on to the quartic form, examine y = ax 4 - bx 2 +(b- a). (We do not worry about missing 
consecutive powers of x because we can translate the above form either right or left and all the con- 
secutive powers of x will appear.) 

Again examine the discriminate knowing that 1 is an automatic root. 

This becomes (- b ) 2 - 4 (a)(b -a) = b 2 - 4 ab + 4a 2 =(b- 2a ) 2 . Thus, the expression ax 4 - bx 2 +(b-a) is 
always factorable regardless of the values for b and 2a. 


Take the derivative, in order to see what values of a and b would make it factorable; this gives 


y = a • x 4 - b • x 2 + (b - a) 


y f = 4a • x 3 - 2b • x 
0 = 2x-(2tf -x 2 -b) 
2x = 0, 2a • x 2 - b = 0 
2 a-x 2 =b 


In order to get integral roots, — must equal a perfect square such as 1, 4, 9, 16. 

2 a 


— = 1,4,9,16,25... 

2 a 

b = 2a, 8a, 18a, 32a, 50a... 


So going back to the original y = ax 4 - bx 2 +(b- a), substitute b = 2a and obtain 
y = ax 4 - 2 ax 2 + a = a(x 4 - 2x 2 + 1) which gives double roots (disregard); letting b = 8a, 
y = ax 4 - 8 ax 2 + la = a{x - 8x 2 + 7) = a{x - 7)(x 2 - 1) which contains irrational roots. 
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Continue the search until rational roots appear (TABLE II). 


TABLE II 


Value of b 

Resulting Polynomial 

Factors of Polynomial 

2a 

ax 4 - 2 ax 2 + a 

a(x 4 - 2x 2 + 1) = a(x 2 - 1) 2 

8a 

ax 4 -8 ax 2 + la 

a(x 4 - 8x 2 + 7) = a(x 2 - l)(x 2 - 7) 

18a 

ax 4 -18 ax 2 +17 a 

a(x 4 -18x 2 + 17) = a(x 2 -l)(x 2 -17) 

32a 

ax 4 -32 ax 2 +31 a 

a(x 4 - 32x 2 +31) = a(x 2 - l)(x 2 - 3 1) 

50a 

ax 4 -50 ax 2 +49 a 

a(x 4 - 50x 2 + 4a) = a(x 2 - l)(x 2 - 49) 


( 5 ) 

Here a • (x 2 - 1) • (x 2 - 49) = a • (x + 1) • (x - 1) • (x - 7) • (x + 7) and the ensuing sequence for TABLE II 

1 , 7, 1 7, 3 1 , 49, . . . 2n 2 - 1 does not yield any perfect squares for relatively small values 
but a new family of quartics is acquired because the derivative is also factorable at n = 7. 

or, 

y - a-(x 4 - 50x 2 + 49) 
y - a- (4x 3 -lOOx) 
y - a- (4x) • (x 2 - 25) 
y - a- (4x) • (x - 5) • (x + 5) 

thus for 

0 = a • (4x) • (x - 5) • (x + 5) 
x = 0, x = 5,x = —5 

TABLE III reflects the new family of quartics and some of the family members. It suffices to say 
that translation left or right can go on indefinitely, but the selection is limited to the list below because 
the numbers become large. 

TABLE III: 

Family 


Translated Form 

x(x + 2)(x - 6)(x + 9) 

(x + l)(x + 3)( x- 5)(x + 9) 
(x + 2)(x + 4)(x - 4)(x +10) 
(x + 3)(x + 5)(x - 3)(x +11) 
(x + 4)(x + 6)(x - 2)(x +12) 
(x + 5)(x + 7)(x - l)(x +13) 
x(x + 6)(x + 8)(x +14) 


Derivative 

4(x + l)(x + 6)(x - 4) 
4(x + 2)(x + 7)(x - 3) 
4(x + 3)(x + 8)(x - 2) 
4(x + 4)(x + 9)(x - 1) 
4x(x + 5)(x +10) 

4(x + 6)(x + 1 l)(x + 1) 
4(x + 7)(x + 12)(x + 2) 


(x - l)(x + l)(x - 7)(x + 7) 
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(7) 

So for integral roots, — - 1 must be a perfect square or 
a 

--1 = 1,4,9,16,25,36,49,64,81,100... 
a 

- = 2,5,10,17,26,37,50,65,82,101... 
a 

or b- 2a,5a,\0a,\7 a,26a,27 a,50a,65a,82a,\0\a... 

For an example: let b = 26 a then 

y - ax 5 - 26ax 3 + 25 ax = ax(x 4 - 26x 2 + 25) = ax(x 2 - 25)(x 2 - 1) = ax(x - 5)(x + 5)(x - l)(x + 1) 
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In terms of the derivative 

y = ax 5 - bx 3 + (b- a)x 
y = 5 ax 4 - 3bx 2 + b-a 

and B 2 - 4 AC = (-3 b) 2 - 4(5 a)(b - a) 

= 9b 2 -20ab + 20a 1 

Finding values a and b for which the above is factorable gives rise to TABLE IV . 


( 8 ) 

TABLE IV 


Value of a 

Resulting Polynomial 

Values of b 

1 

9b 2 -20b + 20 

-1,1,2 

2 

9b 2 - 40b + SO 

-2,1,2,4,11 

3 

9b 2 - 60b + 180 

-3, 3, 6 

4 

9b 2 -806 + 320 

-4,2,4, 8,11,22 

5 

9b 2 -1006 + 500 

-22, -5, 5, 10 

6 

9b 2 -\20b + 120 

-6, 3, 6, 12, 33 

7 

9b 2 -140& + 980 

-19, -7, -2, 7, 11, 14, 29 

8 

9b 2 -1606 + 1280 

4, 8, 16, 22,44 

9 

9b 2 -1806 + 1620 

2, 9, 11, 18, 29 

10 

9b 2 - 200b + 2000 

5, 10, 20, 38 


etc. 


Unfortunately, negative values of b in TABLE IV create complex roots in the original polynomial 
and multiples of “a ” create double or triple roots. Thus our choices are quite limited when stating 
new families. Hence, there are no monic polynomials (a = \) but there are families where a ^ 1... —> . 
A good example here is a = 9 and b = 29 and the analysis gives the following results: 

Here 

y = 9x 5 - 29 x 3 + 20x = x(9x 4 - 29x 2 + 20) = x(9x 4 - 9x 2 - 20x 2 + 20) = x(9x 2 (x 2 - 1) - 20(x 2 - 1)) 



3 


3 
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( 9 ) 


for the derivative 

y = 9x 5 - 29 x 3 + 20x 
y' = 45x 4 - 87x 2 + 20 
y' = 45x 4 - 15x 2 - 12x 2 + 20 
0 = 15x 2 (3x 2 -5)-4(3x 2 -5) 

0 = (3x 2 - 5)(1 5x 2 - 4) 

o = (x41 - 41 )(xV3 + 41 )(xVl5 - 2)(x^ + 2 ) 


V 5 4$ Vf 5 

-4s V 3 — v/l 5 

Xl ~4^'<S~ 3 

41 41~ 3 

2 4l5 2V15 

12 

<N 

1 

<N 

1 

X3 ~4Ts vn _ 15 

4 ~4Ts 4Ts~ 15 


The reader is left to develop other such families where 2, 3, 4 ... are automatic roots. These can 
all be achieved using the remainder theorem. 


Summary; 

The sole objective of this paper was to present both families of polynomials which are 
easily factored and methods by which to determine them. The process is even easier if we 
plant 1, or 2, or 3 as roots making the original function factorable. And, since the only two 
numerical coefficients are a and b, we can determine values of b easily while holding “a ’’ 
constant. Limiting the tables to small numbers is another effort to make successful curve 
sketching of polynomial functions more accessible to students of mathematics. 
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Last summer started an extended email exchange between Heinrich Ludwig and me. 

His mails are given below together with a summary in English. His Snub Cube turned 
out to be the starting point of an exciting journey into the world of polyhedrons. You 
are friendly invited to follow us. 

Von: Heinrich Ludwig [heinrich.ludwig@gym-raubling.de] 

Gesendet: Sonntag, 24. Juli 2011 11:27 

Betreff: DERIVE: Gefullte Darstellung eines Vielecks im 3D-Fenster 
Lieber Herr Bohm, 

meine Frage klingt simpel, ich habe aber in alien DNLs, die in pdf-Version verfugbar sind, keinen 
Hinweis auf eine Antwort gefunden. 

Vielleicht konnten Sie mir weiterhelfen? 

Wie stellt man mit DERIVE 6 ein n-Eck (Eckenzahl n > 4) dessen Ecken dreidimensionale Koordina- 
ten haben, als Flache gefarbt im 3D-Fenster dar? (Naturlich miissen die Eckpunkte in einer Ebene 
liegen.) 

Fur n=3 und n=4 habe ich eine Moglichkeit gefunden. Man bildet aus den Seiten des n-Ecks einen 
Vektor, wobei der Endpunkt der einen Strecke der Anfangspunkt der folgenden Strecke ist. Fur n=5 
klappt das nicht mehr, denn intern behandelt DERIVE dieses Konstrukt so, dass zwei aufeinander 
folgende Strecken als Dreieck aufgefasst werden, welches flachig gefarbt dann dargestellt wird. 

Ich habe auch die Funktion POLY FILL(u) ausprobiert, aber leider wird damit aus jedem Vieleck ein 
’’Zebra”, denn POLY FILL unterteilt das Vieleck in Dreiecke, farbt diese verschieden ein und zeich- 
net auch die Diagonalen des Vielecks, welche es in Dreiecke unterteilen, als schwarze Striche ein. 

Weil der Rechenaufwand zur Priifung, ob ein Vieleck eben ist, betrachtlich sein diirfte (vermutlich n - 
2 Kreuzprodukte und n - 3 Tests auf Kollinearitat der Normalenvektoren) nehme ich nach all meinem 
Ausprobieren an, dass DERIVE nur auf das Einfarben von Dreiecken eingerichtet ist. Ich wollte ar- 
chimedische Korper nicht nur mit ihrem Kantengeriist, sondern als gefullte Korper darstellen. Da ware 
es schon gewesen, wenn man auch 4-, 5-, 6-, 8- und 10-Ecke hatte einfarben konnen. 

Eine Stelle im DNL fand ich, von der ich mir Hilfe versprochen hatte: Im DNL #77 haben Sie die 
Abituraufgabe mit DERIVE gelost, die spater den medienwirksamen Titel ’’Oktaeder des Grauens” 
erhalten hat. Das abgestumpfte Oktaeder haben Sie dann 3 -dimensional dargestellt. 

Ich habe die Datei horror.dfw, die in gezippter Form dem Newsletter 
beilag, mit DERIVE geladen und die entsprechenden Zeilen, die die 
POLYGON FILL-Funktion enthalten, zum Zeichnen verwendet. Und 
jetzt kommt das, was mir unverstandlich ist: Bei Ihnen zeichnet POLY- 
GON FILL schone einfarbige Polygone, bei mir aber diese ’’Zebras”. 

Was lauft da verkehrt? Gibt es eine Einstellung, die DERIVES Farb- 
wechselschema steuert? Ich finde nirgends etwas Geeignetes in den Me- 
nus. Damit Sie sehen, was bei mir passiert, hange ich ein Bild zu meinem 
Ergebnis an dieses Mail. 

Ich hore jetzt auf mit meinen Nachforschungen und wiirde mich recht freuen, wenn Sie mir den ent- 
scheidenden Tipp geben konnten. 

Mit besten GriiBen 
Ihr Heinrich Ludwig 
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Heinrich Ludwig’s problem from 24 July 2011 summarized: 

How can I represent a planar polygon (number of vertices n > 4) in 3D-space as a coloured 
area? 

I found a way for n = 3 and n = 4: 

One has to form a vector consisting of the edges of the n-gon with the endpoint of one side 
being the start point of the next one. But this does not work for n > 4, because DERIVE inter- 
prets two subsequent segments as a triangle. 

Then I tried the function POLY_FILL(u) but unfortunately every polygon appears as a “ze- 
bra”. POLY_FILL divides the polygon into triangles, dyes them in different colours and plots 
all sides of all triangles in black. 

I assume that DERIVE needs some calculation for testing if the points are lying in a common 
plane etc (which needs some cross products) and that the function is only intended for col- 
ouring triangles. I wanted to represent Archimedean solids not only as a wireframe model but 
as a filled solid. It would be great to dye 4-, 5-, 6- 8- and 10-gons. 

I found one article in a DNL which promised some support. 

In DNL#77 you presented a truncated octahedron (“Octahedron of Horror”) using function 
POLYGON_FILL. I tried this function and I cannot understand what happened then: Your 
polygons are pretty one-coloured - and I received again “zebras”. What am I doing wrong? Is 
there a DERIVE setting which controls the colour change. 

I will stop now my investigations and would be happy if you could provide any hint. 

Sincerely Yours 
Heinrich Ludwig 


See my answer 
Dear Mr. Ludwig, 

Your “zebra” will change to a “grey donkey” or a “brown stallion” if you switch off option 
“Automatically change color of new plots” under Options > Display > Color. Plotting the bor- 
ders needs some “manual work”: 

Perform a right mouse click on each triangle and the Edit-popup window will open: Choose 
Scheme: Custom, then choose the colour for the polygon and set the same colour for the 
mesh lines. 

Performing this procedure for all triangles should lead to the desired result. 

Best regards 
Josef 
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Von: Heinrich Ludwig [mailto:heinrich.ludwig@gym-raubling.de] 

Gesendet: Dienstag, 26. Juli 2011 23:15 
Betreff: kleines Dankeschon fur Ihre Hilfe 

Lieber Herr Bohm, 

als kleines Dankeschon fur Ihre freundliche Hilfestellung schicke ich Ihnen, wo fur ich das Einfarben 
von Flachen verwenden wollte: einen archimedischen Korper. Mein Liebling ist der Cubus Simus. Er 
tanzt so schon aus der Reihe, ist eben nicht zentralsymmetrisch wie fast alle anderen. AuBerdem zer- 
fallen die Dreiecke unter den Seitenflachen des Cubus Simus topologisch in zwei Klassen, auch eine 
Besonderheit, die nur noch ein weiterer archimedischer Korper (das Dodekahedron Simum) hat. 

In den DNL-Heften, sowohl in den alten, die Sie neu als pdf aufgelegt haben, wie auch in den neueren, 
habe ich nichts Explizites zu archimedischen Korpern gefunden, insbesondere keine Koordinatensatze 
fur die Ecken. Ich kann mir nicht vorstellen, dass es keinen schonen Beitrag dazu gegeben hat. Er wird 
wohl in einem der Hefte 26 bis 52 stehen, fur die es keine pdf-Version gibt und die ich darum nicht 
kenne. Weil so ein Beitrag sehr wahrscheinlich schon erschienen ist, arbeite ich die Sache nicht weiter 
aus. Fur Schuler ist der Cubus Simus eine zu harte Nuss. Nicht einmal meinem Schuler, der sich in 
seiner Seminararbeit mit der 3D-Darstellung verschiedener Korper beschaftigt, traue ich das zu. (Es 
gibt ja geniigend viele einfachere Korper auch noch.) Aber weil das Wetter in den letzten Tagen gar so 
bedauerlich war, habe ich selber ein Objekt zum Spielen gebraucht. 

Der erste Schritt zur Konstruktion des Korpers war es, die Koordinaten der Ecken anzugeben. Ich 
wusste nur soviel iiber sie, dass sie auf die Seitenflachen eines Wiirfels gelegt werden konnen und 
dazu kannte ich die Symmetrien des Cubus Simus. Mit Hilfe zweier Parameter k und m, deren Werte 
noch zu bestimmen waren, konnte ich die Koordinaten formulieren. Der zweite Schritt bestand dann 
darin, die Ecken richtig miteinander zu verbinden. Dazu muss der Graph des Polyeders aufgestellt 
werden. Hat man beides, dann ist es nur noch eine FleiBaufgabe, die 38 Seitenflachen einzutippen. 

Da ich's ohnehin nur fur mich aufbewahre, geniigt mir eine handschriftliche Skizze, wie man zu den 
Eckpunktskoordinaten kommen kann. Den Graphen habe ich mit Herrn Hohenwarters Geogebra ge- 
zeichnet (siehe Seite 25, Josef). 

Es wird an unserer Schule gem verwendet und ich selber nutze es auch recht gerne. 

Herzliche GriiBe 
Heinrich Ludwig 

Heinrich Ludwig’s mail from 26 July 2011 summarized: 

As a little expression of thank I am sending an Archimedean solid. This was what I needed 
colouring areas for. My favourite solid is the Cubus Simus. It is pretty extraordinary because 
it is not central symmetric like almost all others. Additionally the faces are divided topologi- 
cally in two classes. This is also a special feature only shared with one other Archimedean 
solid, the Dodecahedron Simum. 

I didn’t find any contribution on Archimedean Solids in earlier DNLs, in particular no sets of 
coordinates for the vertices. I wonder if there has no respective article appeared up to now. 
The Cubus Simus is a hard nut for students. Not even one of my students who works on a 
paper containing 3D-representations of various solids would be able to do this. (Fortunately 
enough there are many other easier solids, too.) 


Because of bad weather in the last days I needed an object to play with. 
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The first step for constructing the solid was to 
express the coordinates of its vertices. All what 
I knew was that they are located on the faces of 
a cube - and I knew the symmetries of the Cu- 
bus Simus. Using two parameters k and m 
(their values had to be found) I could express 
the coordinates, (see my sketch). 

The second step was to connect the vertices 
(24) in the right way (60 edges forming 6 
squares and 32 triangles). This needs compos- 
ing the graph of the polyhedron (see the Geo- 
Gebra graph). The rest is some diligent work: 
entering the 38 faces of the solid. 




The GeoGebra Graph 


Der kronrandige Wurfel, das abgeschragte Hexaeder, 

Cubus Simus, Snub Cube 

Heinrich Ludwig, Germany 

Der Korper ist dem Wurfel mit den Ecken [±1, ±1, ±1] einbeschrieben. Alle Ecken des Korpers liegen 
auf den Seitenflachen des Wiirfels. Die Koordinaten der 24 Ecken A bis X sind ±1, ±k oder ±m. Die 
genauen Werte fur k und m werden unten hergeleitet. Es ist auch interessant zu sehen, wie sich der 
Korper verandert, wenn man k und m mit Schiebem variiert. Man kann auch ein Kubooktaeder oder 
ein kleines Rhombikubooktaeder erzeugen. 

The solid is inscribed a cube with vertices [±1 , ±1 , ±1], All 24 vertices - A through X - lie in 
the faces of the cube with coordinates ±1 , ±k or ±m. The exact values for k and m are de- 
rived below. It is interesting to observe how the solid changes varying k and m by applying 
sliders. It is possible to create a cuboctahedron or a small rhombicuboctahedron. 
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Cuboctahedron 


Small Rhombicuboctahedron 


Here are the 24 vertices: 


#6: 

[A 

:= [1, 

k , hi ] , B 

:= [1, -k, -ill], C 

:= [1, -in. 

k], D := 

[1, in, -k]] 


#7: 

[ E 

:= [HI, 

1, k], F 

:= [-I", 1, -k], G 

:= [k, 1, 

-111], H := 

[-k, 1, m]] 


#8: 

[I 

:= [k, 

m , 1], 3 

:= [ k , -in, 1], K 

:= [-ill, k, 

1], L := 

[in, -k, 1]] 


#9: 

[M 

:= [-1, 

-k, ni], 

N := [-1, 

k, -ni], 

0 := [-1, 

in, k], P 

:= [-1, -hi, 

-k]] 

#10: 

[Q 

:= [-l», 

-1, k], 

R := [ill, 

-1, -k], 

S := [k, - 

-1, in], T 

:= [-k , -1, 

-in]] 

#11: 

[U 

:= [-k, 

in, -l], 

-XL 
1 1 

> 

-m, -l], 

W :: [-111, 

-k, -1], 

X :z [ill, k, 

-i]] 


These are the 6 squares followed by 32 triangles: 



A 

D 


E 

H 


I 

1 



D 

B 


H 

F 


1 

J 


ADBC :z 

B 

C 

, EHFG := 

F 

G 

1 1 

1 — 1 

_l 

1— 1 

J 

K 


- 

C 

A 


G 

E 


K 

I 

- 

- 

H 

P 


' Q 

T 


U 

X 

- 


P 

N 


T 

R 


X 

V 


HPNO := 

N 

0 

, QTRS := 

R 

S 

, UXVW := 

V 

w 



0 

M 


S 

Q . 


w 

u 




' A 

E ' 


' C 

A ' 


' C 

I ' 


' I 

E ' 


AEI := 

E 

I 

, CAI := 

A 

I 

, CIL := 

I 

L 

, IEK := 

E 

K 


- 

. I 

A . 


. I 

C . 


. L 

C . 


. K 

I . 

- 

- 

" E 

H ' 


' A 

G ‘ 


‘ D 

G ' 


" L 

C ' 

- 

EHK := 

H 

K 

, AGE := 

G 

E 

, DGA := 

G 

A 

, LCS := 

C 

S 



. K 

E . 


. E 

A . 


. A 

D . 


. S 

L . 
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' D Q ] [ K 0 ] D G 

JQM := Q M , KOH := 0 H , DGX := G X , RVB := 

H J J H K J X D 

' P W ] [ R W ] [ R T ] [ P U 

PWT := W T , RWV := W V , RTW := T W , PUW :: U W 

. T P 1 [vrJ [wrJ [WP 

' P N ] [XUl [UN] [FN 

PNU := N U , XUF := U F , UNF := N F , FNH := N H 

. U P J [ F X J [fuJ [HF 

H N ] [ J K ] [DO] Q H 

HNO := N 0 , JKO := K 0 , DOM := 0 M , QHT := H T 

„ 0 H J |_ 0 J J M D J |_ T Q 

T H ] [VD] [VX] G F 

THP := H P , VDB := D B , VXD := X D , GFX := F X 

. P T 1 IbvJ [dvJ [XG 

L D ] [QS] [SB] [ R B 

LDQ := D Q , QSL := S L , SBC := B C , RBS := B S 

.qlJ [lqJ [csj S R 

Calculation of appropriate values for k and nr. 

The edges CA, Al and 1C must be of equal length. 

2 2 2 
| C — A | = 2 • ( k + m ) 

2 2 2 
I A — 1 1 =2 •(k - k*(m + 1) + m - m + 1) 

2 2 2 
1 1 - C | = 2 • C k - 2-k + 2-m + 1) 

2 2 2 2 
SOLVE( | A - I| = |C - A| a | I — C | = |C - A| , [k, m], Real) 
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2 2 2 2 
SOLVE( | A - I| = |C - A| a | I — C | = |C - A| , [k, m], Real) 


k = 


1/3 1/3 1/3 

2- (233 - 39-733) + 2- (39-733 + 233) + 2- (6-733 - 26) - 2- 

18 


1/3 1/3 1/3 

(6-733 +26) - 6 (6-733 - 26) (6-733 + 26) 

a m = - + 


1 

3 

k = 0.5436890126 a m = 0.2955977425 

Plotting all squares and triangles gives a wonderful representation of the snub cube. 



Hint: You can highlight and plot all expressions from ADBC:= to RBS:= in one single step. 

Comment: The values for k and m and closely related with the Tribonacci Constant c, which 
is the real solution of equation x 3 — x * 2 - x - 1 = 0. (c « 1 .839287) 

(Additional Comment: The Tribonacci Sequence is defined as a n = a n - 1 + a n ~2 + a„_ 3 for n > 4 
with a-i = a 2 = 1 , a 3 = 2. 

1, 1,2,4, 7, 13, 24,44, ... 

The ratio of adjacent elements tends to the Tribonacci Constant and the numbers of the se- 
quence increase asymptotically to a n ~ c n . 

(See also my additional comments at the end of this contribution, Josef) 
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3 2 

SOLVE (X -X - X - 1 = 0, X) 
1/3 


1/3 


(19 - 3-733) (3-733 + 19) 1 

X — — — + + L « 

6 6 3 


(27-711 + 57-73) 


1/3 


(57-73 - 27-711) 


1/3 


1/3 1/3 

(19 - 3-733) (3-733 + 19) 1 

v X = - - + + 

6 6 3 


1/3 1/3 ^ 

(57-73 - 27-711) (27-711 + 57-73) 


(19 - 3-733) 


1/3 


v X = 


1/3 

(3-733 + 19) 1 

+ 

3 3 


X = -0.4196433776 + 0 . 6062907292 - L v X = -0.4196433776 - 0 .6062907292 - L v 


1.839286755 


Von: Josef Bohm [nojo.boehm@pgv.at] 

Gesendet: Samstag, 30. Juli 2011 13:25 

Dear Mr. Ludwig, 

I couldn’t resist colouring your pretty Cubus Simus, too. 

I did it in another way: First of all I produced the wireframe model of the solid (using your 
data together with GeoGebra-Graph). Then I copied and pasted the figure from the 3D-Plot 
Window into the Algebra Window. I converted this picture to a Paintbrush graph (Edit > Con- 
vert Picture Object > Paintbrush Picture). Then I could fill all the faces in Paintbrush. The 
edges remain visible. As you can see I chose the visibility from the other side. 



JB 
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Von: Heinrich Ludwig [mailto:heinrich.ludwig@gym-raubling.de] 

Gesendet: Dienstag, 9. August 2011 13:11 

Lieber Herr Bohm, 

es freut mich, dass meine Spielereien Sie angeregt haben, damit weiter zu experimentieren. Ich habe 
auch schon daran gedacht, das Einfarben einem Programm zu iiberlassen, das sich auf graphische Be- 
arbeitung besser versteht als Derive. Liefert Derive das Kantenmodell des Korpers, dann hat ein Gra- 
phikprogramm alle Information, um die weitere Gestaltung zu ubemehmen. Naturlich verliert man 
dabei, dass man auf einfache Weise den Augenpunkt wandem lassen kann, aber fur eine statische Ab- 
bildung des Korpers ist das nicht wichtig. 

Den von Ihnen beschriebenen Weg kann man noch ein wenig abkiirzen, indem man den Graphikinhalt 
des 3D-Fensters mit Strg+C unmittelbar in die Windows-Zwischenablage kopiert und deren Inhalt 
dann vom Graphikprogramm einliest. Als solches bevorzuge ich GIMP. Es ist das am meisten ausge- 
reifte OpenSource-Produkt auf dem Markt. 

Wussten Sie schon, dass 

f 1/4 

2 .(75 - 1 ) 

ATAN 

( 11 • 7T ^ f 29*7T ^ C f B • 7T ^ 

8. COS -COS .y - 4. 72 . SIN - 7(5 - 75 ) + 3-72 

l l 120 J l 120 J L l 20 J 

Fur Derive war die Vereinfachung des Terms leider ein zu harter Brocken. Wie kommt man aber auf 
so einen verriickten Term? Mich hat die Frage beschaftigt (das Wetter war allzu traurig in letzter Zeit), 
wie gut archimedische Korper rollen konnen. Das ist sicher eine eher physikalische Frage. Mein ma- 
thematischer Ansatz war, den Raumwinkel an der Ecke des Polyeders zu betrachten. Nirgends konnte 
ich dazu eine Auflistung der Werte fur die 5+13 archimedischen Korper finden, also habe ich Derive 
anpacken lassen. Derive musste die Formel von L'Huilier auswerten. Und dann kam fur das Hexakisi- 
kosaeder (4,6,10) das o.g. Ungetiim heraus. Jetzt ware es freilich interessant, was die groBen Compu- 
teralgebrasysteme (Mathematica, Maple, etc.) leisten. 

Schaffen diese die Vereinfachung? Und was leistet der Derive-Nachfolger Nspire? 

Herzliche GriiBe 
Heinrich Ludwig 



Heinrich Ludwig’s mail from 9 August 2011 summarized: 

Dear Mr. Bohm 

The procedure which you described can be shortened by copying the contents of the 3D-plot 
window into the clipboard (Ctrl+C) and then pasting into a graphics program. I prefer GIMP. 
For me it seems to be the best available open source program on the market. 



B • 7T 
8 
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It was a work too hard for DERIVE to simplify. But does one encounter such a crazy expres- 
sion? I was busy treating the problem (weather was too bad in the last time) how good 
Archimedean solids can roll along. I know that this is more a physics question. My mathe- 
matical approach was to investigate the space angle in an edge of the solid. As I could not 
find any list of these angles for the 5 + 1 3 Archimedean solids I engaged DERIVE to do the 
job. DERIVE had to evaluate the formulae of L’Huilier. And then there was the Hexakis Ico- 
sahedron. It is the source of the onster from above. I would be very interested what the „big“ 
CAS like MATHEMATICA and Maple are able to do with it. 

Can they perform the simplification? And what about Tl-Nspire? 


Dear Mr. Ludwig, 

Many thanks for your extended letter. Don’t underestimate DERIVE. It needs only setting the 
right mode for treating trig expressions: 


ATAN 


1/4 

2 -C7B - 1) 


8- COS 


11*7T 


• COS 


v 120 y 
Trigonometry := Collect 
( 7(72 + 2 )^| 

ASIN 


29- 7T 


120 


■7 


4-72 -SIN 


3*7T 


20 


- 7(5 - 75 ) + 3.72 


It is also interesting simplifying this expression piecewise: 


11-7T ^1 


8- COS 


120 


29- tt } 


• COS 


120 


= 4 . COS 


7 


4 . cos 


3-tt 


20 


^2 


+ 2 


- 4. 72. SIN 


3 *7T 


20 


3 *7T 


+ 2 


20 


\\ 


- 7(5 - 75) + 3.72 




3/4 1/4 3/4 1/4 

75.2 - 50 - 2 +2 

1/4 

2 -(75 - 1 ) 

3/4 1/4 3/4 1/4 

75.2 - 50 - 2 +2 


= 72+1 


ATAN (72 + 1) = 


8 
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or: 


r 

- 4.72 -SIN 

" 3 «TT ' 

\ 

- 7(5 - 75) + 3.72 • 

r 

■ 4* COS 

r 3 *7T ' 

\ 

+ 2 

V 

. 20 , 

J 

\ 

, 20 > 

J 


- 6-710 + 8-75 + 18- 72 - 24 
1/4 

2 -C75 - 1) 

= 72+1 

7(- 6.710 + 8.75 + 18- 72 - 24) 


You could also perform “stepwise simplification” in order to get some insight how DERIVE 
works internally. 


I add picture of the Hexakis Icosahedron. It 
is one of the Catalan solids which are dual to 
the Archimedean solids. 

The H.l. is the dual of the Great Rhombico- 
sidodecahedron. 

(Source: the interesting kauernet-website 
from below.) 


Some additional comments: 

(1 ) Inserting sliders for k and m provides more insight in the influence of the parameters. 







Challenge 1: Which combination of k and m create the solids given on page 26? 

(2) Another definition of the snub cube: 

The solid whose vertices are the 24 points on the surface of a sphere for which the 
smallest distance between any two is as great as possible. 
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(3) Heinrich Ludwig had in the head line of his DERIVE file the expression Pentagonal 
Icositetrahedron. This is the “dual” solid to the snub cube. What are dual polyhedra? 

The dual of a polyhedron can be found by connecting the midpoints of the edges sur- 
rounding each vertex to a polygon and then constructing its tangent polygon (of the cir- 
cumcircle). The numbers of vertices and faces interchange. The P.l. has 24 faces and 
38 vertices. 

The dual solids of the Archimedean solids are the Catalan solids. 

Challenge 2: Find the DERIVE model of the Pentagonal Icositetrahedron. 

(4) For our German readers: a snub nose is a “Stupsnaschen”. 

For all readers: “Simus” is a latin word for “abgeplattet” (German) and “flattened, oblate” 

(English). 

(5) The Tribonacci-sequence 

trib_c(20, 1, 1, 2) = [1, 1, 2, 4, 7, 13, 24, 44, 81, 149, 274, 504, 

927, 1705, 3136, 5768, 10609, 19513, 35890, 66012] 

(trib_c(100, 1, 1, 2)) = 98079530178586034536500564 

100 

1/3 1/3 

(19 - 3-^33) (3-733 + 19) 1 

c := + + 

3 3 3 


SSI. 839286755214161132551852564653286600424 


(trib_c(21, 1, 1, 2)) 

21 

(trib_c(20, 1, 1, 2)) 

20 

(trib_c(101, 1, 1, 2)) 

101 

(trib_c(100, 1, 1, 2)) 

100 


1.839286796340059383142458946858147003575 


1.839286755214161132551852564653286600424 


The next page illustrates the properties of the Tribonacci constant. Finding a function 
for generating the Tribonacci sequence is left for the reader. 
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r (trib_c(x + 1, 1, 1, 2)) 


TABLE 


x + 1 


(trib_c(x 5 i, i, 2)) 


x, 1, 20 


c 

2 ■ ■ 


!■ 


5 

x 

c 

TABLE ((trib_c(x s i, i, 2)) s x, 1, 50) 

x 



(6) A (small) selection of (many) interesting websites. (I didn’t add the Wikipedia sites 
which are also highly recommended.) 

http : //arxi v. org/pdf/ 1202.0172vl .pdf 

http://www.ewetel.net/-neubauer/raumgeometrie/dateien/basteln/cubus%20simus%20basteln.pdf 
http://mathworld.wolfram.com/PentagonalIcositetrahedron.html 
http://kauernet.de/hendes/html/h306 catal.htm 


10 15 20 


x 

LN(c ) 

TABLE C(LN(trib_c(x s 1, 1, 2))) , x 5 1, 100) 



http://mathworld.wolfram.com/LHuiliersTheorem.html 

http://www.calcuhis.conLkonfuciy.asp?tda=dt&t=966&fs=l%27huilier%27s+theorem 

(7) Finally, many thanks to Heinrich Ludwig for his inspiring contribution. 

(8) Most of my “knowledge” and some pictures presented above are from the internet. 

Josef 
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Statistics 

Visualising and Simulating Dynamically with Tl-Nspire 3.1 

Guido Herweyers 
KHBO Campus Oostende 
quido.herwevers@khbo.be 


Introduction 

“Measuring is Knowing” is more important than ever nowadays. Wednesday 20 October 2010 
was the first “World Statistics Day” initiated by the UN. At this occasion the NIS (Nationaal 
Instituut voor Statistiek) [1] 2 published a special extended version of the core numbers 2009 
with Belgium placed into an European perspective 171 . 

Technology proceeds and offers new possibilities for teaching. This is for sure for statistics. 
Actual data can be visualised quickly in various ways, simulations of random samples and 
experiments can be investigated. Doing so statistical concepts can be introduced early placing 
the data into the foreground and the formulae in the background. The intention of this paper is 
to demonstrate a possible way offering concrete examples. The statistical basic ideas are not 
defined here. The necessary background can be found in numerous statistics textbooks. 

Part 1 of this paper is to introduce the possibilities of Tl-Nspire (version 3 . 1 ) [3] , containing 
mainly descriptive statistics (2 nd and 3 rd grade in Belgium with two to four lessons mathemat- 
ics a week). We are working with data from the web, too. 

Part 2 deals with simulating with Tl-Nspire. How to select a simple random sample from a 
population? How does the normal distribution appear as probability distribution for the sam- 
ple mean observing graphically the variability of more than 1000 sample means? 

Part 3 shows how to discover new probability distributions by simulation, starting with the 
uniform distribution on the interval [0,1]. We also treat the t-distribution. 

Part 4 demonstrates a possible approach for hypothesis testing supported by simulations. 

This paper can be downloaded (in its Flemish original) from the T Vlaanderen website, 

www.t3vlaanderen.be , 

together with the Tl-Nspire files and all the data appearing in this paper. 


[1] NIS (Nationaal Instituut voor Statistiek) is the biggest statistical organisation of the Belgian 
government: http://statbel.feov.be/nl 

[2] Kemcijfers 2009, Belgie in een Europees perspectief 
http://economie.feov.be/nl/modules/pressrelease/statistieken/generale/world statistics day.jsp 


[3] Information about Tl-Nspire: www.education.ti.com 
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Part 1 : Descriptive Statistics 

(1) Qualitative ungrouped data 

Open TI-Nspire, close the welcome screen and insert a page with the Lists & Spreadsheet Applica- 
tion (which is one of seven applications offered): 



Example 1: 

Given are the eye colours of 15 persons: green, blue, blue, grey, blue, green, brown, blue, grey, brown, 
blue, green, blue, blue and grey. 

— » Name the first column A as eye_colour. 

—> Enter the colours (begin in cell A1 ) one after another - under quotes because the data are strings. 



Document 1 - Tl-Nspire™ CAS Student Software 

File Edit View Insert Tools Window Help 



fS In: 

sert - ^ 


1 

► 


Documents ToolboM 

II 

o 

-C: 

Q 

nr 

1 i 8 M rf 


E 


1- 

Ml 

£= 

Oj 

E 

3 

o 

O 

Q 





Lists & Spreadsheet £ 

^ 1:Actions ► 

§g2:lnsert ► 









— 

1 ' 3 : 3: Data ► ' 

1 :Generate Sequence 
2\ Data Capture ► 

3:FiU 

4 : C te a r Data 

^Summary Plot 

X 4: Statistics ► J 

5:T able ► 



1 

3:Quick Graph 


n — ' 




11 
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Select the whole column (left mouse click on letter A), then choose the first menu (Docu- 

ment Tools) in the Document Toolbox. Choose option 3: Data > Quick Graph. The window will 
split and you will receive a Dot Plot in the right window. 



A eye_colour 

B 

■ 

■ 

♦ 





1 

green 




2 

blue 




3 

blue 




4 

grey 




5 

blue 




6 

green 




7 

brown 




8 

blue 




9 

grey 




10 

brown 




11 

blue 




12 

green 




13 

blue 




14 

blue 




15 

grey 



— a 

A 

eye_colour 

◄ ► 


1:Bar Chart 


3:Sort 


blue 


brown green 
eye_colour == 


grey 


— » Choose the right window (left mouse click on the window) and change the diagram to a Pie Chart. 
Another (right) mouse click on the pie offers presenting the percentages (Show all Labels). Hovering 
the mouse over the sectors of the pie opens a small window showing details. 



A eye_colour 

■ 

c 

♦ 

i 




green 




2 

3 

4 

blue 

blue 

grey 



5 

blue 



6 

green 



7 

brown 



8 

blue 


l 

9 

grey 


10 

11 

brown 

blue 


A1 "green" 



green 


bro., 



blue 


eye_colour 


The context menu (right mouse click) enables quick performing appropriate actions. You can close a 
window by selecting it (left mouse click), pressing Ctrl + K followed by pressing the Del(ete)-key. 
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— » Change the diagram to a Bar Chart (right mouse click in the right window). 

Shift the category brown to the right (left mouse click on the label brown and drag it to the right). 



3 (20.<*») 3 (20.0°o) 


2C13.3H) 

I 


green grey brown 

eye colour 


Right mouse click pop up menu offers three 
ways of sorting the bars. 


1 :Dot Chart 
2:Pie Chart 


3:Hide All Labels 

|4:Sort ► 

1 : List Order 
2:Value Order 
3:Alpha Order 



— » Add another column containing the gender of the persons, and then split the category of the eye 
colour with respect to the second category gender. (Right click on the X- variable eye_colour, select 
option 4: Split Categories per Variable and select gender.) 


hi 

A eye_colour ^ gender 

C 

D a 

♦ 




1 

green boy 



2 

3 

4 

5 

blue girl 

blue girl 

grey girl 

blue boy 



6 

green girl 



7 

brown boy 



8 

blue boy 



9 

grey boy 



10 

11 

12 

brown boy 

blue girl 

green girl 



13 

blue boy 



14 

blue girl 



15 

grey girl 


— H 

A1 "green" 



4- 


& 
c 

Q) 

D 
cr 

,?2H 


0 J 


rsi 

S' 

m 



blue 


green 


grey 


o 

d 

(N 

d 


brown 


eye_colour / gender 
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—> Insert a Calculator page via Insert from the Menu bar. 


Documents Toolbox 

> ' Problem 


j ' Page Strg+I 



Data & Statistics 


Graphs 

Geometry 

Lists & Spreadsheet 
Jj Data & Statistics 
Notes 

Vernier DataGuest 

[Si. 1_:Plot Type 
; 2: PLot Properties 
3: Actio ns 
4: Analyze 
tHg 5:Window/Zoom 

Image... 

Video... 

PublishView™ Text Box 


|0 Program Editor ► 

(0 Sensor Console Strg+D 


-^►Type eye_colour followed by pressing the ENTER-key. The variable eye_colour is a list of 
strings. The command countif(eye_colour,?=,”blue”) gives the number of appearances of “blue” in 
this list. 


eye_colour 

{ " green" , " blue" , " blue" , 11 grey 11 , 11 blue" , "green" , " brown " , " blue " , 
c0untIf|^e_c0/0Mr ? ?= ,, blue ,, j 7 


Hint: The size of the characters can be adapted by using the slider bar on the bottom of the 


screen: 


Document View: 


Boldness: 


L±_ 


=t:t= 


..i.j 


(2) Qualitative grouped data 

Example 2: 

On the NIS website one can find the distribution of the use of ground in Belgium (2009) 


Bodemgebruik Belgie 2009 (in km 2 ) 
bron: NIS 

Landbouw 

10387,34 

Permanente weiden en grasland 

4963,39 

Bossen 

6970,57 

Andere en water 

2156,42 

Bebouwde percelen 

6050,19 


(agriculture, permanent pastures and grassland, woods, others and water, tilled areas) 


— » Insert a new spreadsheet page and fill in the data. 
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Insert a new spreadsheet page and fill in the data. 


Documents Toolbox 

| 1 9 SB 


Lists & Spreadsheet 


E+3 IjActions 
^ 2: Insert 


3:Data 


X 4: Statistics 
5:Table 


km c 


87.3 


~a ~aQ 


1 :Generate Sequence 
2: Data Capture > 

3: FiU 

4: Clear Data 


5:Summary Plot 


6: Quick Graph 



A ground 

B square_km 

c 

D 

■ 

♦ 







Agriculture 


10387.3 


2 

Grassland 

4963.39 


3 

4 

Woods 

Others & Water 

6970.57 

2156.42 


5 

Tilled Area 

6050.19 



X List ground 


Summary List: square_km ▼ 


Display On: Split Page ▼ 


OK Cancel 



A ground 

B square_km 

C 

♦ 




1 

Agriculture 

10387.3 


2 

Grassland 

4963.39 


3 

4 

Woods 

Others & Water 

6970.57 

2156.42 


5 

Tilled Area 

6050.19 


6 

7 



8 




9 




10 




11 

12 




1 -> 




m 





A5 "Tilled Area" 

* 

► 


10000 


8000 


l 6000 




Grassland 
4963.39 ( 16 . 3 %) 


ground 


—> Activate the right window. Bring the bars to a decreasing order by the menu Actions > Sort > 
Value Order. You can also shift the bars respectively. 
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Then change the bar chart to a pie chart via the Plot Type menu (or by a right mouse click on one of 
the bars). 


Data & Statistics 

* 

[Sf 1_: Plot Type 

► 

Wl 2: Plot Properties 

► 

| 3:Actions ► 

% 4:Analyze 

► 

5:Window/Zoom 

► 


|re_km 

C E 
" 1 



)387.3 


3Q 



fcyd 1: Remove 
2: Hide Text 

Abl 3:lnsertText 
4:lnsert Slider 
Djj 5: Select All Points 

Vs 6: Select Image 


10000 - 

8000- 

;i 6000- 

i 

r 4000 
2000 



-p 

O' 

O 

'Sf 

+ 



The statistics diagrams appear in a Data & Statistics Application. You could also start by inserting 
this application into a separate page, then use the Plot Properties menu, select Add X Variable 
(ground) and finally select Add Y Summary List (square_km). 


Insert t 




E Problem 
( 1 Page 


Strg+I 


Calculator 

Graphs 

Geometry 

a Lists & Spreadsheet 


ffu Data & Statistics 


Notes 

J Vernier DataGuest 
|e| Image... 

jr Video... 

PublishView™ Text Box 
Program Editor 


in Sensor Console 


Strg+D 


Data & Statistics 


1_: Plot Type 


: 2: Plot Properties 


3: Actio ns 
" n "- u - 4: Analyze 
^ 5:Window/Zoom 


1 :Connect Data Points 
E4 2:Histogram Properties ► 

3:Extend Box Plot Whiskers 
4:Show All Labels 


6:Add X Summary List 
7:Remove X Variable 
’ 8:Add Y Variable 
9:Add Y Summary List 
A: Remove Y Variable 


B:Force Numeric X 
C:Force Numeric Y 
D: Split Categories by Variable 
E:Remove Split Variable 


F:Clear All 


But it works really faster by a rmcl (right mouse click) on Click to add variable on the window bot- 
tom and then choosing Add X Variable with Summary List. 


Click to add 

: i 

1 :Add Variable 



2:Add X Variable with Summary List || 
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Save the file as statistics (Menu File > Save as). 


Open menu Page Sorter 
as Qualitative Variables by a rmcl on Problem 1. 


in the Documents Toolbox. Rename the standard name Problem 1 


, statistics - Tl-Nspire™ CAS Student Software 


@3 New Tl-Nspire Document Strg+N 

@ New PublishView 1 " Document Strg*Umschalt*N 

Open Document... Strg+O 

Close Strg+W 


B Save Document 


Save to Handheld.. 
Convert to 
Export 


Print... 


Settings 


View Copyright Information 


Recent Documents 


Strg+S 


Strg+P 


Documents Toolbox 


T Problem 


Cut 

Strg+X 

Copy 

Strg+C 

Paste 

Strg+V 

f ' Page 

Strg+I 

® Problem 


Rename F2 

If Delete 

Entf I 


Documents Toolbox 

it HI e " tf 


T Qualitative Variables 


[* 



One file can contain various problems, with each problem containing more pages (which can be se- 


lected in the Page Sorter 


l. Each page can consist of maximum four applications (seven applica- 


tions are available). Each application has its own menu which is available under 
ments Toolbox. 



in the Docu- 


Within one problem a variable has the same value on all pages. Variable x in one problem is inde- 
pendent of variable x in another problem. 


Example 3: 

Belgium consists of three districts (constituent states): the Flemish district (Vlaams Gewest), the Brus- 
sels Capital Town district (Brussels Hoofstedelijk Gewest) and the Wallonian district (Waals district). 
The German speaking minority (Duitstalige gemeenschap) is part of the Waals district. 

The population numbers from 2005 to 2010 were (according NIS): 


Gewest 

2005 

2006 

2007 

2008 

2009 

2010 

Brussels H. Gew. 

1.006.749 

1.018.804 

1.031.215 

1.048.491 

1.048.491 

1.089.538 

Vlaams Gewest 

6.043.161 

6.078.600 

6.117.440 

6.161.600 

6.208.877 

6.251.983 

Waals Gewest 
waarvan 
Duitstalige gem. 

3.395.942 

3.413.978 

3.435.879 

3.456.775 

3.475.671 

3.498.384 

75.512 

73.119 

73.675 

74.169 

74.540 

75.222 


—> Create a bar chart comparing district populations of the years 2005 and 2010. 
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Hint: you can add different Add Y Summary Lists with one X-variable given (via Plot Type in the 


Data & Statistics Documents Toolbox 



A gewest 

B ly2005 

®y2006 

D y2007 

E ly2008 

F ly2009 

G y2010 

♦ 








1 

Brussels 

1006749 

1018804 

1031215 

1 04849lJ 

1068532J 

1089538 

2 

Vlaams 

6043161 

6078600 

6117440 

6161600 

6208877 

6251983 

3 

Waals 

3395942 

3413978 

3435879 

3456775 

3475671 

3498384[ 



Our USER Forum 

On 07/01/12 13:53, Giuseppe Omaghi wrote: 

I am using Derive 6.10 for Windows and I would like to know if it is possible to plot the absolute 
value of the zeta function in a range. 

Thank you very much. 

Giuseppe 

Hi Giuseppe, 

try ABS(ZETA(x+i*y)) as a 3-d plot. You will get a surface in R A 3 defined by 
[x,y,ABS(ZETA(x+i*y))]. I hope this is what you were looking for. 


Stefan Welke 
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With ABS(ZETA(x+i*y)) I get the error message:Too many variables.With z=ABS(ZETA(0.5+i*y)) I 
get a surface in 3D. I would like to plot ABS(ZETA(0.5+i*y)) on a xy plane (x-axis corresponds to the 
values of y and y-axis to the absolute value of ZETA(0.5+i*y)). 

Is it possible? 

Thank you very much. 

Giuseppe 


Hello Giuseppe, 

did you use the "i" from the menu in the bottom? If not, then "i" from the keyboard will be interpreted 
as a variable, and then there are three, which is one too much. Alternatively try ABS(ZETA(x+#i*y)) 
as input. This will work correctly. I assume, that you will now see a surface different from that of your 
previous attempt. 

Cordiali saluti, Stefan Welke 



I see a surface different from previous examples. Thank you. 
And for my attempt to plot a graph in xy-plane? 

Thank you very much. 

Saluti dalfltalia. 

Giuseppe 

|£(0. 5 + l-x) 



(symmetric wrt y-axis) 
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TIME 2012 - The Keynotes, Lectures and Workshops 

List of presentations (in alphabetical order) of accepted proposals. It is very 
likely that more will follow. All continents will be present in Estonia in July. 

A Stepwise Graphical Approach for Teaching Automated Theorem Proving in Engi- 
neering 

A Tool for Evaluation of Rationality of Algebraic Transformations 
Adapting the Assessment to the Instrumental Genesis 

Adopting Tl-Nspire CAS Technology (Handheld and Software), a Campus-Wide Ex- 
perience 

An Evaluation of Students: Experiences in Technology-based Courses in Mathemat- 
ics and Statistics 

Attracted by (Strange) Attractors 

Being a Math Teacher - Looking through the Blog Posts of Student Teachers (Key- 
note) 

CAS in the Classroom: Yesterday, Today, and Tomorrow (Keynote) 

CAS or not CAS? that is the question when using technology in math education 

Changing Assessment Methods: New Rules, New Roles 

Compound Damped Pendulum: an Energy Comparison 

Computer Algebra as an Educational Tool 

Critical Issues of Effective CAS Utilization (Keynote) 

Deriving Big Formulas with Derive and What Happened Then 

Determining the Proof Skills of 2nd Year Preservice Mathematics Teachers by 
Means of GeoGebra 

Development of Pedagogical Tools for Teaching Mathematics from the PC Lab to 
Smartphone (Keynote) 

Dynamic Systems - From Simulation Software to CAS (WS) 

Dynamic Systems - Many Competences - Presented Dynamically 
Experiences of 15 Years Using Graphic Calculators in Saxonia 
Finding Bioinformatics Data and DERIVE Programs to Process it (WS) 

Find your very own "attractive" Attractor 

Finnish Upper Secondary Students' Collaborative Processes in Learning Statistics in 
CSCL Environment 

From Inductive Reasoning to Proof by Induction with Geometry Expressions (WS) 

From Vienna to Tartu: A 10 Years Tour using DERIVE in TIME (Keynote) 

Hidden Markov Models and an Introduction to Their Uses in Bioinformatics 

High School Students' Acquisition of Knowledge and Skills through Self-Organization 
and Collaboration 

Integration of Digital Tools into the Mathematics Classroom: A Challenge for Prepar- 
ing and Supporting the Teacher 
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Invariant Region Method (IRM) Applied to Two Dimensional Predator-Prey Systems 
Learning Mathematics through the Surrounding World 

Making Learning Math Addictive - Localization of the Famous Open Source Website 
Khan Academy in Estonia 

Mathematical Induction and CAS 

Mathematical Modelling with Geometry Expressions and Tl-Nspire CAS 
Math eTextbooks 

Media Diversity Project - E-Learning Paths in Maths Education 
Modeling the Probability of a Binary Outcome 

Motion Sensors in Mathematics Teaching: Learning Tools for Understanding General 
Math Concepts? 

Notice Mathematics around You 

Programming-like Activities for Calculus and Differential Equations Classes 

Prospective Teachers' Curriculum Strategies in CAS-based Written Lessons 

Representations of Solids and Surfaces within the Tl-Nspire Environment 

Systems of Computer Algebra and Dynamic Geometry as Tools of Mathematical 
Investigation 

Teaching conics through mathematical modelling (WS) 

Technology in Maths (a Queensland Perspective) 

The Effects of Using a Calculator on Reducing Mathematics Anxiety 

Types of Mistakes and Dynamics of the Occurrence of Mistakes Solving Linear 
Equations Using Computer Program T-algebra 

Using Derive To Generate DNA Fractal Representations 

Using Dynamic Documents to Enhance Student Understanding 

Using Fourier Series to Analyse Mass Imperfections in Vibratory Gyroscopes 

Using Tl-Nspire CAS Technology in Teaching Engineering Mathematics: ODE's 

Using Tl-Nspire CAS Technology in Teaching Engineering Mathematics: Calculus 

Visualizing and Understanding Statistical Concepts with Tl-Nspire 

What was Special in DERIVE? 

When do Computer Algebra Systems Offer Unexpected Answers to School Equa- 
tions? 

Why Trying to Avoid Using Complex Numbers? 

DUG-Meeting 2012 

For more Information and Registration: 
http://www.time2012.ut.ee/ 


